Precise stellarator quasi-symmetry can be achieved with electromagnetic coils

Magnetic fields with quasi-symmetry are known to provide good confinement of charged particles and plasmas, but the extent to which quasi-symmetry can be achieved in practice has remained an open question. Recent work [M. Landreman and E. Paul, Phys. Rev. Lett. 128, 035001, 2022] reports the discovery of toroidal magnetic fields that are quasi-symmetric to orders-of-magnitude higher precision than previously known fields. We show that these fields can be accurately produced using electromagnetic coils of only moderate engineering complexity, that is, coils that have low curvature and that are sufficiently separated from each other. Our results demonstrate that these new quasi-symmetric fields are relevant for applications requiring the confinement of energetic charged particles for long time scales, such as nuclear fusion. The coils’ length plays an important role for how well the quasi-symmetric fields can be approximated. For the longest coil set considered and a mean field strength of 1 T, the departure from quasi-symmetry is of the order of Earth’s magnetic field. Additionally, we find that magnetic surfaces extend far outside the plasma boundary used by Landreman and Paul, providing confinement far from the core. Simulations confirm that the magnetic fields generated by the new coils confine particles with high kinetic energy substantially longer than previously known coil configurations. In particular, when scaled to a reactor, the best found configuration loses only 0.04% of energetic particles born at midradius when following guiding center trajectories for 200 ms.

Controlled nuclear fusion is a promising candidate to satisfy rising electricity needs while avoiding carbon emission into the atmosphere. To generate electricity from fusion, a plasma needs to be maintained at extremely high temperatures over long time scales, which requires excellent particle confinement. This is typically achieved using powerful magnets. Tokamaks rely on a toroidally axisymmetric system of magnetic coils to achieve good confinement, at the cost of requiring a plasma current to generate a significant fraction of the magnetic field. This plasma current is challenging to drive in steady-state operation (1), and can be the source of disruptive instabilities (2,3). In contrast, the nonaxisymmetric coil systems of stellarators can generate a confining magnetic field in the absence of plasma currents, relieving many of the challenges for continuous, disruptionfree operation. However, the lack of axisymmetry implies that neither nested magnetic flux surfaces nor particle confinement is guaranteed (4). This motivates the need for a generalization of axisymmetry to the stellarator context, called quasi-symmetry.
A magnetic field is said to satisfy quasi-symmetry if there exists an invariant direction for the field strength B = |B| in a certain coordinate system (5). This condition leads to the conservation of canonical angular momentum, which, in turn, implies the remarkable property that such fields are guaranteed to confine charged particles without requiring plasma currents. However, it remains an open question whether three-dimensional magnetic fields that are perfectly quasi-symmetric over a volumetric region exist. Recently, using numerical optimization, Landreman and Paul (6) found vacuum magnetic fields that satisfy the quasi-symmetry property in toroidal geometries to a very high precision. Magnetic fields that are quasi-symmetric to such a high degree have not been discovered before. Simulations confirm excellent confinement properties even for particles with a large kinetic energy. Naturally, the question arises whether these magnetic fields can be accurately produced by a set of practical electromagnets, which would be a first step toward a new generation of fusion experiments. In this report, we show that such magnets, in fact, exist. Specifically, we find coils producing fields whose deviation from perfect quasisymmetry is more than four orders of magnitude smaller than the mean field strength. For a 1-T mean field, comparable to many stellarator experiments, this error can be on the order of Earth's magnetic field. Simulations confirm particle confinement comparable to the ideal fields discovered in ref. 6. This shifts the discovery of these highly quasi-symmetric fields from one of theoretical interest to one of practical relevance for fusion experiments and other confinement applications (7). . We note that the rectangular cross-section of the coils is for visualization only; in all computations for this manuscript, the coils are each approximated by a single wire of infinitesimal thickness. We conducted numerical experiments in which we approximated finite-thickness coils using multiple filaments, and obtained very similar results, not shown here.

Results and Discussion
We solve a constrained optimization problem to obtain coils that approximate the two quasi-axisymmetric fields of ref. 6, which we refer to as "QA-LP" and "QA+Well-LP," reflecting the absence or presence of a magnetic well. The design space consists of four distinct modular coils, which results in 16 coils in total after applying symmetries. The length of coils impacts the quality of the magnetic field approximation strongly, and hence we compute coil sets of different coil length L max and compare their performance, for example, "QA+Well[20]" refers to the coil set which approximates the "QA+Well-LP" configuration and for which the four modular coils have combined length 20 m. For both "QA-LP" and "QA+Well-LP," as we allow longer coils, the field induced by the coils becomes a better approximation of the target field. Values of L max are relative to an average major radius of 1 m. Fig. 1 shows coils obtained by choosing L max = 18 m and L max = 24 m in the approximation to QA-LP, the relative normal magnetic field B · n/|B| on the surface S, and a Poincaré plot for L max = 24 m. Here B is computed using the Biot-Savart law, and, if the normal magnetic field is exactly zero, then the field matches that discovered in ref. 6 everywhere in the volume contained by the surface. For the shortest coils, we observe the largest normal magnetic field on the surface, reaching values of up to 3.3 × 10 −3 . Its oscillatory nature is caused by the discrete nature of the electromagnetic coils and their close proximity to the surface. Longer coils enable a larger distance between the surface and the magnets, thereby reducing these discrete effects, and more accurately reproducing the target magnetic field, with a relative normal magnetic field of, at most, 1.6 × 10 −4 for the QA[24] configuration.
The Poincaré plot demonstrates the existence of nested magnetic surfaces in the entire volume. This feature is necessary (but not sufficient) for confinement, since charged particle trajectories are tangent to these surfaces in the limit of low energy.
Quasi-axisymmetric fields are characterized by the following property: When parametrized using Boozer coordinates φ, θ, the field strength on each magnetic surface only depends on the angle θ (5,8). When this property is satisfied exactly, even highly energetic collisionless particles are confined over long time scales, and collisional transport is minimized. In Fig. 2C We quantify this statement by performing a Fourier transformation of |B| and then studying the magnitude of those Fourier coefficients that break quasi-symmetry; that is, we write |B(s, θ, φ)| = m,n B m,n (s) cos(mθ − nφ) and then consider terms with n = 0. Here s indexes each surface by the normalized toroidal magnetic flux it encloses, and we plot the largest symmetry-breaking Fourier mode on each surface in Fig. 2A. As coil length is increased, the symmetry-breaking error is reduced significantly, and approximates those of the fields discovered in ref. 6. In fact, we are able to achieve errors smaller than Earth's magnetic field for the longest set of coils when the mean field strength is 1 T. For comparison, gray curves show the symmetry-breaking amplitudes for the eight previous quasi-symmetric configurations in figure 1 of ref. 6.
Finally, scaling the configurations to the mean field and minor radius of the ARIES-CS reactor (9), we compute collisionless guiding center trajectories for alpha particles initialized with a kinetic energy of 3.5 MeV, as by-products of a deuterium-tritium fusion reaction, on the surface with normalized toroidal flux s = 0.25 (half radius). Particle losses are shown in Fig. 2B. For L max = 24 m, the performance is nearly indistinguishable from the target equilibrium, and, for the QA+Well[24] configuration, less than 0.04% of particles are lost after 0.2 s, a typical time for the alphas to thermalize with the main plasma. This coil set is shown in Fig. 2D. Performance is only slightly worse for L max = 22 m, but is poor for the coils with L max = 18 m. For comparison, gray curves show calculations for the nine previous stellarator configurations of figure 6A in the SI Appendix of ref. 6, similarly scaled. We note that the thermal collisional transport magnitude <6 × 10 −8 for QA [24]. These values are orders of magnitude below the values for other optimized stellarators [e.g., ∼10 −3 for Wendelstein 7-X (4)], and are so small that collisional fluxes would be negligible compared to turbulent transport. This is the standard situation for tokamaks, but is unusual and desirable for stellarators, since collisional and turbulent losses are additive (4).
In conclusion, we have shown that the magnetic fields of ref. 6 can be produced very accurately using coils, making these fields practically relevant for stellarators. As a result, exceptionally good confinement of particle trajectories and remarkably small thermal collisional transport can be achieved. While longer coils are required for optimal performance, these coils are not particularly complex as measured in terms of curvature and coil-to-coil separation.

Materials and Methods
The electromagnetic coils were optimized using the SIMSOPT software (11). Similar to the approach in ref. 12, each coil is modeled as a closed, smooth curve in R 3 and represented using a Fourier series, truncated at order 16. Given a magnetic surface S (obtained from ref. 6), the objective that we minimize is given by where B is the field induced by the coils. If f B = 0, then the induced field is exactly equal to the target field up to a scaling factor.
Finding coils that minimize this objective is an ill-posed problem, so we require additional regularization. In practice, it is desirable to have coils that are not too long, avoid high curvature, and are well separated. In this work, we enforce constraints on the curvature (κmax ≤ 5 m −1 ), the mean squared curvature (κmsc ≤ 5 m −2 ), and the distance between coils (d min ≥ 0.1 m), and we vary the constraint on the total length of the magnetic coils (Lmax ∈ {18, 20, 22, and 24 m}). The units quoted in the above assume a major radius scaled to 1 m. We refer to SI Appendix for more detail on the exact implementation of these constraints.
The optimization objective has 399 parameters for the coils and is highly nonconvex. The minimization uses the L-BFGS-B algorithm with analytic gradients. To remedy the possibility of the optimizer being stuck in a local minimum, we start the optimization from eight different initial coil sets and choose the best minimizer as measured by the objective. The code is highly optimized and parallelized and is publicly available at https://github.com/ florianwechsung/CoilsForPreciseQS. Using eight cores of an Intel Xeon Platinum 8268 CPU, solving an individual optimization problem takes ∼30 min. Data Availability. Code and optimization results have been deposited in GitHub (https://github.com/florianwechsung/CoilsForPreciseQS) and Zenodo (10.5281/zenodo.5975323) (13).